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We study to what extent the measured elliptic flow at RHIC constrains viscous deviations 
from ideal hydrodynamics. We solve a toy model where only transverse momenta are 
thermalized while the system undergoes longitudinal free-streaming. We show that RHIC 
data exclude such a model and thus require fast 3- dimensional thermalization. 

1. MOTIVATION 

One of the most important discoveries at RHIC so far has been the strong elliptic flow 
generated in non-central collisions [ 0, 0, |^ . For not too peripheral collisions the measured 
elliptic flow coefficient V2{p±) at midrapidity almost exhausts the hydrodynamic limit [ 
^ ^ ^, up to transverse momenta p± <2 GeV, and in this entire domain the dependence 
of f 2 on the mass of the emitted particles [ ||, accurately follows the hydrodynamically 
predicted pattern [ ^. This means that the distribution of the momenta of well over 99% 
of the emitted particles is accurately described by (ideal) hydrodynamics. 

Why is this so important? The initial transverse momentum distribution of the particles 
generated by the colliding nuclei is locally isotropic. Only their spatial distribution in the 
transverse plane is initially deformed (for b^O). Interactions among the produced quanta 
are required to transfer this spatial anisotropy to momentum space. A non-vanishing f 2 is 
thus an unambiguous signature for reinteractions in the produced matter, and the observed 
large V2 values prove that the initial parton liberation process is separated from the 
experimentally observed final state by a violently interacting stage of dynamical evolution. 
Consequently, there is a priori little reason to expect that calculable properties of the 
early matter formed in the reaction zone immediately after nuclear impact (such as, 
for example, the central rapidity density and the shapes of the rapidity and transverse 
momentum distributions of the produced gluons [ ^) have any direct relationship with 
the corresponding experimentally observed values. The evolution of spectral shapes due 
to the strong rescattering must be taken into account, and even though adiabatic cooling 
and the build-up of collective flow move the spectral slopes in opposite directions it is 
highly unlikely that these effects cancel completely. Elliptic flow itself is an example for 
a qualitative change of the spectra between particle formation and decoupling. 

Microscopic studies [ |ll| show that V2{pi_) is a monotonic function of the mean 
free path, X^{ap)~^; the hydrodynamic limit is approached from below as A— s>0. The 
observation that at RHIC V2{p±) almost exhausts the hydrodynamic limit thus appears 
to force the conclusion [ |l], ^ |, 0, |12| that the fireballs formed in Au+Au collisions at 
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RHIC thermalize fast and efficiently. Since the creation of elliptic flow is driven by the 
spatial deformation of the reaction zone which quickly decreases either spontaneously via 
free-streaming [ |^ or (more rapidly) as a result of the flow anisotropy itself [ |T3| , f 2 is 
sensitive to the very early collision stage. Within hydrodynamics the RHIC V2 data can 
only be reproduced by assuming thermalization at times rth < 1 fm/ c (we use rth = 0.6 fm/c 
[ 1^, 1^, 0). As the hydrodynamic limit of V2 can only be reached from below and the 
data almost exhaust it, we can use the hydrodynamic model to estimate the energy 
density at thermalization. One finds values [ |^ at least an order of magnitude above 
the critical one for hadronization, e ^ 1 GeV/fm^. The thermalized state formed ~ 1 fm/c 
after nuclear impact must therefore have been a quark-gluon plasma which, according to 
the hydrodynamical simulations, lives for several fm/c before hadrons ffist appear. Unless 
the presented chain of arguments leading to the conclusion of early thermalization can be 
broken, the implication is unavoidable that at RHIC a well-developed quark-gluon plasma 
has been created. 

The purpose of the work reported here [ |14| is to try to poke holes into the early 
thermalization argument. It was already pointed out by Ollitrault [ ^ that V2 is sensitive 
to the stiffness of the equation of state (EOS) of the thermalized matter and increases 
monotonically with the sound velocity c^ = dP/ de. Might it be possible to trade off 
thermalization against a stiffer EOS? The concrete idea studied by us [ |14[ was that 
perhaps the thermalization of the (on average much larger) longitudinal momenta of 
the liberated partons takes longer than transverse thermalization, resulting in a smaller 
longitudinal than transverse thermal pressure. As shown by Teaney at this conference, this 
is similar to the expected effects from shear viscosity on the longitudinal hydrodynamic 



expansion [|T5|. If we exaggerate a bit and idealize the model by assuming collisionless free- 
streaming with boost-invariant initial conditions in the longitudinal direction coupled with 
complete local thermalization of the transverse momenta, we can write down an analytic 
expression for the phase-space distribution function in terms of macroscopic parameters 
for which we can analytically derive macroscopic equations of motion which generalize the 
usual ideal hydrodynamic equations. In this idealization there is no longitudinal pressure 
at all; all the hydrodynamic work goes in the transverse direction. Can we generate larger 
elliptic flow in this way? If yes, the RHIC data would no longer saturate the theoretical 
limit and the fast thermalization argument would break down. We'll see that it doesn't. 



2. A TRANSVERSALLY THERMALIZED MODEL (TTHM) 

With boost-invariant initial conditions, assuming that all particles originate at t = 
from z = (i.e. infinitely Lorentz contracted colliding nuclei), an appropriate ansatz for 
the phase-space distribution of massless gluons is 

/(x,fc,t) = ^5(y-r,)-l-^,^. (1) 

Here r = y/t'^—z'^, 77 = | In j^, and y = ^\n fn^f^. We use longitudinal boost-invariance to 
parametrize the local flow as = '-f±{cosh.7], v±, sinh?]). The transverse flow velocity v± 

with 7_|_ = l/-y/l— f j_ and the "transverse temperature" T are functions of the longitudinal 
proper time r and the transverse position x±, but independent of space-time rapidity rj. 
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Due to the free- streaming constraint 6{y—ri) the exponent of the Bose distribution reduces 
to k-u = 'y±k±{l — k±-v_i_), showing thermahzation of only the transverse momenta. 

Inserting this ansatz into the kinetic definition of the gluon energy momentum tensor, 
T^'^{x, t) = Ug J f{x, k, t) {vg = 16 is the gluon spin-color degeneracy), one finds 

T^'^ = (e + P±)u^'u'' - P± Ig^"" + n'^n" + (w'^m'^+m^M'^)] (2) 



with the equation of state (EOS) e{x±,T) = 2P±{x±,t) = z/y(ro/T)(7r^T^(a;_L, r)/60). 
The difference between longitudinal and transverse momenta requires two additional vec- 
tors to decompose T^^", = (sinh?7, 0, 0, coshr^) and m'^ = 7_l(w_l cosh 77, v±, v± sinh?]) with 
n-m = n-u = m-u = 0. The resulting extra terms relative to the ideal fluid decomposition 



can be viewed as viscous corrections [ [T5[. The EOS is consistent with a vanishing trace 
of T^'^ for massless particles and the absence of longitudinal pressure; the corresponding 
sound velocity Cg = 1/a/2 is larger than for an ideal gluon gas, i.e. the EOS is stiffer. 
With the decomposition (|^) the macroscopic equations of motion d^T^'^ = become 

+ — + ViT^^'' = 0, z/ = 0,l,2,3, J = 1,2,3. (3) 

OT T 

The only difference from the analogous ideal hydrodynamic equations [ ^ is the absence 
of a term —P/t on the right hand side of the z/ = equation, reflecting work done by 
the longitudinal pressure. Standard hydrodynamic codes for the transverse expansion of 
systems with longitudinal boost invariance can thus be used to solve these equations. 



3. RESULTS 

Lacking longitudinal pressure, the energy density in TTHM decreases more slowly with 
r than in ideal hydrodynamics (HDM). However, the larger transverse pressure performs 
more work in the transverse directions, resulting in stronger radial and elliptic flow. As a 
consequence, initial conditions which in HDM lead to a consistent description of the data 
produce in TTHM much too fiat transverse momentum spectra. The reproduction of the 
central collision data via TTHM thus requires retuned initial conditions. This retuning 
is facilitated by the numerical observation [ |l^ that for massless particles the TTHM 
dynamics leads to completely time-independent transverse momentum spectra; transverse 
flow buildup thus exactly compensates for cooling! We have no analytical explanation for 
this fact. It implies that the final slope of the transverse momentum spectrum is given by 
the initial temperature; steeper final spectra thus require a reduced initial temperature. 
To maintain the measured normalization dN/dy one then must also increase the initial 
longitudinal volume, by increasing the starting time tq for the TTHM dynamics. This 
implies late transverse thermahzation (and, of course, an even later longitudinal one). 

To obtain the same final spectra and radial flow in TTHM as in HDM tuned to RHIC 
data, we find that we should start the transverse expansion about 10 times later with 
a ~ 15 — 20 times lower initial energy density. The total remaining time until freeze- 
out is thereby shortened considerably, and even though we do not take into account the 
decrease of the spatial deformation of the reaction zone by transverse free-streaming prior 
to the onset of transverse TTHM dynamics, we find that now considerably less elliptic 
flow (only about half as much as in HDM) is generated (Fig. 1). As the data almost 
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exhaust the HDM values, TTHM underpredicts them by about a factor two, both for the 
Px-integrated elhptic flow f 2 (Fig. la) and for the p^-slope of the differential elliptic flow 
V2{p±) (Fig- lb). The TTHM model is thus experimentally excluded; only a model such 
as HDM with a considerable degree of early longitudinal thermalization can describe the 
RHIC measurements. 
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Figure 1. V2 as a function of impact parameter h (left) and of for two values of h (right) 
from TTHM (solid) and HDM (dashed). For details see Ref. [ 0]. 

In closing we note that for massless bosons f2(p_L— »0) approaches a non-zero value 
(Fig. lb). This reflects the singular nature of the Bose distribution at p_L = [ We 
suggest to use this as a test for thermalization of small-p^ gluons in parton cascades and 
of direct photons in experiment. 
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